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Abstract 

Employing Hilbert-Schmidt measure, we explicitly compute and analyze a number of determi- 
nantal product (bivariate) moments |p| fc | / cr PT | n , k, n = 0, 1, 2, 3, . . ., PT denoting partial transpose, 
for both generic (9-dimensional) two-rebit (a = |) and generic (15-dimensional) two-qubit (a = 1) 
density matrices p. The results are, then, incorporated by Dunkl into a general formula (Ap- 
pendix |D6j), parameterized by k,n and a, with the case a = 2, presumptively corresponding to 
generic (27-dimensional) quaternionic systems. Holding the Dyson-index-like parameter a fixed, 
the induced univariate moments l/O^* 7 " | and \p PT \ n are inputted into a Legendre-polynomial- 
based (least-squares) probability-distribution reconstruction algorithm of Provost (Mathematica 
J., 9, 727 (2005)), yielding a-specific separability probability estimates. Since, as the number of 
inputted moments grows, estimates based on the variable |/?||p PT | strongly decrease, while ones 
employing \p PT \ strongly increase (and converge faster), the gaps between upper and lower esti- 
mates diminish, yielding sharper and sharper bounds. Remarkably, for a = 2, with the use of 
2,325 moments, a separability-probability lower-bound 0.999999987 as large as k. 0.0804954 is 
found. For a = 1, based on 2,415 moments, a lower bound results that is 0.999997066 times as 
large as J| ~ 0.242424, a (simpler still) fractional value that had previously been conjectured (J. 
Phys. A, 40, 14279 (2007)). Furthermore, for a = ^, employing 3,310 moments, the lower bound 
is 0.999955 times as large as || = 0.453125, a rational value previously considered (J. Phys. A, 
43, 195302 (2010)). 
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I. INTRODUCTION 



In a much cited paper pQ, Zyczkowski, Horodecki, Sanpera and Lewenstein expanded 
upon "three main reasons" -"philosophical", "practical" and "physical" -for attempting to 
evaluate the probability that mixed states of composite quantum systems are separable 
in nature. Pursuing such a research agenda, it was conjectured [21 sec. IX]-based on 
"a confluence of numerical and theoretical results" -that the separability probabilities of 
generic (15-dimensional) two-qubit and (9-dimensional) two-rebit quantum systems, in terms 
of the Hilbert-Schmidt /Euclidean/fiat (HS) measures [31 03, are ^ « 0.242424 and § « 
0.470588, respectively. In this study, we shall avail ourselves of newly-proposed formulas 
of Dunkl (Appendix D) for (bivariate) moments of products of determinants of density 
matrices (p) and of their partial transposes (p PT ) [3 E] to investigate these hypotheses from 
a novel perspective, as well as extend our analyses beyond the strictly two-rebit and two- 
qubit frameworks. (To be fully explicit, we note here that both [symmetric] two-rebit and 
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[Hermitian] two-qubit 4x4 density matrices p have unit trace and nonnegative eigenvalues, 
while their partial transposes p PT can be obtained by transposing in place the four 2x2 
blocks of p. The Hilbert-Schmidt metric-from which the corresponding measure can, of 
course, be derived-is defined by the line element squared, ^Tr[(dp) 2 ] pfl, eq. (14.29)].) 

Reconstructions of probability distributions based on these product moment formulas of 



Dunkl do prove to be highly supportive of the specific HS two-qubit conjecture (sec. VIII B), 



while definitively ruling out its two-rebit counterpart (sec. VIII A), but emphatically not a 
later advanced value of || = 0.453125 [7J p. 6]. Extending these analyses from the real 
(a = |) and complex (a = 1) cases to the (presumptively, since we lack relevant computer- 
algebraic determinantal moment calculations) generic (27-dimensional) quaternionic (a = 2) 
instance [8HTT] , in which the off-diagonal entries of the 4x4 density matrices can be quater- 
nions, we find that the value || « 0.0804954 fits our moment-based computations, may 



we say, amazingly well (sec. VIII E). Nevertheless, the apparently formidable challenges of 
rigorously proving the determinantal moment formulas of Dunkl and/or the conjectured 
simple fractional separability probabilities certainly remain. (To again be explicit, the only 
rigorously demonstrated results reported in this paper are those we have been able to ob- 
tain through computer algebraic [Mathematica] methods-using the Cholesky-decomposition 
parameterization of p-for the moments of |p| fc |p PT | n for n = 1,2, . . . , 13 for the two-rebit 
systems and n — 1, 2, 3, 4 for the two-qubit systems [sec. [Ilj, and n = 1 for their qubit-qutrit 
[6 x 6] counterparts [sec. VI , as well as n = 1, . . . , 10 for minimally degenerate two-rebit 
systems [sec. VII] . Aside from the presentation and discussion of these results, the paper is 



concerned with the [unproven] generalization to arbitrary n by Dunkl of these specific results, 
and its apparent successful application in probability-distribution reconstruction procedures 
[sec. VIII . This latter step is taken in order to examine anew and extend certain conjectures 



as to the specific values of the separability probabilities, the properties of which were first 
investigated by Zyczkowski, Horodecki, Sanpera and Lewenstein pQ.) 

In marked contrast to the finite-dimensional focus in this study on 2 x 2 quantum systems 
(and, marginally, on 2x3 systems [sec. |VI| ), let us note the (asymptotically-based) conclusion 
of Ye that "the probability of finding separable quantum states within quantum states is 
extremely small and the Peres-Horodecki PPT criterion as tools to detect separability is 
imprecise for large N, in the sense of both Hilbert-Schmidt and Bures volumes" [121 p. 14]. 
(The Bures distance measures the length of a curve within the cone of positive operators 
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on the Hilbert space [U sec. 9.4], while the Bures volume of the set of mixed states is 
remarkably equal to the volume of an (iV 2 — l)-dimensional hypersphere of radius | [H p. 
351].) Also, contrastingly, to the predominantly "nondegenerate/full-rank" objectives here 



(cf. sec. VII), Ruskai and Werner have demonstrated that 'bipartite states of low rank are 



almost surely entangled" [T3] . 

II. DENSITY-MATRIX DETERMINANTAL PRODUCT MOMENTS 

Let us begin our investigation into the indicated statistical aspects of the "geometry of 
quantum states" jU [H] by noting the two following special cases-which will be extended 
in certain bivariate directions-of the (univariate determinantal moment) formulas [T5] [eq. 
(3.2)] (cf. [El Theorem 4]): 

/U|*\ _ 01 ./>2 fc r(2fc + 2)r(2fc + 4h m 

\\P\ )2-reMt/HS- M5 { A T(4k + 10) J W 

and 

m, », m - 108972864000 r <* + ML + f < * + 3 > F <* + 4 >, (2) 

/2-qubit/HS r(4(fc + 4)) 

k = 0, 1,2, . . . The bracket notation () is employed to denote expected value, while p in- 
dicates a generic (symmetric) two-rebit or generic (Hermitian) two-qubit (4 x 4) density 
matrix. The expectation is taken with respect to the probability distribution determined 
by the Hilbert-Schmidt/Euclidean/flat metric on either the 9-dimensional space of generic 
two-rebit or 15-dimensional space of generic two-qubit systems [3111]. 

At the outset of our study, we were able to compute seventeen (thirteen two-rebit and 
four two-qubit) non-trivial (bivariate) extensions of these two formulas, involving now in 
addition to \p\, the quantum-theoretically important determinant |p PT |. (The nonnegativity 
of |p PT |-as a corollary of the celebrated Peres- Horodeccy results [5j|6]-constitutes a necessary 
and sufficient condition for separability/disentanglement, when p is a 4 x 4 density matrix 
[ITJ HE] ) At this point of our presentation, we note that three of these seventeen extensions 
are expressible-incorporating as the last factors on their right-hand sides, the two formulas 
above (Q, @)-as 

, fe| PT,x (fc-l)(fc(2fc + ll) + 16) , fcx 

PWP \h-reUt/HS 32(A; + 3)(4fc + ll)(4A; + 13) KlPl ^-rebU/HS> W 
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, k PT 2 , _ k(k(k(k(Ak(k + 12) + 203) + 368) + 709) + 2940) + 4860 , k , 

\\P\ \P I /2-rebit/HS ~ 102 4(A; + 3)(fc + 4)(4Jfe + ll)(4Jfe + 13) (Ak + 15) (Ak + 17) ' ' Pl ^-rebit/HS 

(4) 

and 

/\n\ k \n PT \\ fc(fc(fc + 6)-l)~42 , fc . 

^ Vir-vMt/HS %(2k + V)(Ak + \7)(Ak + IV) Km )2 -i uhit / HS ~ K ) 

These three new formulas were, initially, established by "brute force" computation-that 
is calculating the first (k = 0, 1, 2, . . . , 15 or so) instances of them, then employing the 
Mathematica command FindSequenceFunction, and verifying the formulas generated on 
still higher values of k. 

Let us note here the ranges of the two variables of central interest, \p\ G [0, ^g] and 
\p PT \ g [— jg, 255]. For various analytical and conventional purposes, it is often convenient 
to have variables defined over the unit interval [0,1]. If we so (linearly) transform the two 
determinantal variables, then the rational factors on the right-hand sides of ^ and Q get 
replaced, respectively, by 

8(fc(£;(34£; + 297) + 867) + 842) 



17(Jfe + 3)(4Jfe + ll)(4Jb + 13) 

and 



(6) 



6A(k(k(k(k(68k(17k + 348) + 200835) + 904492) + 2279781) + 3048904) + 1689900) 

289(£; + 3)(fc + 4)(4fc + ll)(4fc + 13)(4fc + 15)(4fc + 17) ' ^ ' 

III. THE MIXED /BALANCED VARIABLE \p\\p PT \ = \pp PT \ 

As a special case (k = 1) of formula (|3]), we obtain the rather remarkable moment 
result, zero, already reported in [19]. The immediate interpretation of this finding is that 
for the generic two-rebit systems, the two determinants \p\ and |p PT | comprise a pair of 
nine-dimensional orthogonal polynomials [2"0Tf2"2"] with respect to Hilbert-Schmidt measure. 
(C. Dunkl has kindly pointed out that orthogonality here does not imply zero correlation. 
The analogous quantity for generic twoqubit systems is not zero, however, but — 457 g 264 .) 
In addition to this first (k = 1) HS zero-moment of the product variable |p||p PT | in the 
two-rebit case, we had been able to compute its higher-order moments, k = 2, . . . , 6. (The 
result for k = 2, that is 569634 g 244g00 , can be obtained by direct application of formula Q.) 
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A. Range of Variable 



The feasible range of the (mixed/balanced) variable is |p||p PT | G [— 110 1 592 , T^J-the lower 
bound of which— 110 1 592 = — 2 _12 3~ 3 . This lower bound, determined by analyzing a general 
convex combination of a Bell state and the fully-mixed state, can be achieved with the 
entangled two-rebit density matrix 



P 



eV2 6V2 
6V2 3 12 

(-1 - Vs) 

(-1 - v/3) § 
V^(-l + v / 3) —£m ^ \ ) 



(8) 



6V2 6V2 

The determinant of p here is ^ (2a/3 — 3) ~ 0.000805732 and that of its partial transpose, 
gig (-3 - 2^) « -0.0112224 (their product being « -9.04225 • 10~ 6 ). Both p and 



p PT here have three identical eigenvalues (A (3 — \/3) ~ 0.105662 for p and ^ (3 + \/3) ~ 
0.394338 for p PT ). The isolated eigenvalues for p and p PT are | (l + \/3) « 0.683013, and 
I (l — \/3) ~ —0.183013, respectively. The purity (index of coincidence jH p. 56]) of (8) 
equals |, so the participation ratio is 2. Its concurrence is | (\/3 — l) ~ 0.366025, while its 
entanglement of formation is [U sec. 15.7] 

„ r , log(2) (log (84 + 48^3)-V31og(3)) 

-Ecompiexfp] = " 7^ v ' ~ ~ 1.21665. (9) 

4 log (1 + 73) log(l + V3) 

(Supportively, Dunkl has noted that the computed zeros in his Gaussian quadrature analyses 
D 5 of the two-rebit case fit well into the known ranges of |p| and |p PT |.) Alternatively, 



sec. 



taking into account the real nature of the entries of p, in the sense of the "foil" theory 
of Caves, Fuchs and Rungta [23], one has a concurrence of [2U eq. (4)] and an 

entanglement of formation [24~l eq. (2)] of 

log(2) (log(1728) + 2v / 6tanh~ 1 (Jfj) 

E real [p} = ^- — — \\ » 6.56825. (10) 

m 61og(6-2x/6)log(2(3 + v / 6)) 

IV. CONTOUR PLOTS OF BIVARIATE PROBABILITY DISTRIBUTIONS 

For the further edification of the reader, we present in Fig. [T] a numerically-generated 
contour plot of the joint Hilbert-Schmidt (bivariate) probability distribution of |p| and |p PT | 
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FIG. 1: Contour plot of the joint Hilbert-Schmidt probability distribution of \p\ (horizontal axis) 
and \p PT \ in the two-rebit case. Larger values appear lighter. The variable ranges are \p\ G [0, ^gg] 



and \p 



PT\ 



j_ j_ 

16' 256 



. One billion random density matrices were employed. 



in the two-rebit case, and in Fig. [2j its two-qubit analogue. (A colorized grayscale output is 
employed, in which larger values appear lighter.) In Fig. [3]is displayed the difference obtained 
by subtracting the second (two-qubit) distribution from the first (two-rebit) distribution. 
(The black curves in all three contour plots appear to be attempts by Mathematica to 
establish the nonzero-zero probability boundaries-which, it would, of course, be of interest 
to explicitly determine/parameterize, if possible-of the joint domain of \p\ and |p PT |.) 

These last three figures are based on Hibert- Schmidt sampling (utilizing Ginibre ensem- 
bles [15]) of random density matrices, using 10, 000 = 100 2 bins. In regard to the two-qubit 
plot, K. Zyzckowski informally wrote: "A high peak in the upper corner means that: a) a 
majority of the entangled states is 'little entangled' (small det(p T )) or rather, they are 'close' 
to the boundary of the set, so one eigenvalue is close to zero, and the determinant is small; 
b) as det(p) is also small, it means that these entangled states live close to the boundary of 
the set of all states (at least one eigenvalue is very small), but this is very much consistent 
with the observation that the center of the convex body of the 2-qubit states is separable 
(so entangled states have to live 'close' to the boundary). Similar reasoning has to hold in 
the real case as well." 
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FIG. 2: Contour plot of the joint Hilbert-Schmidt probability distribution of \p\ (horizontal axis) 
and \p PT \ in the two-qubit case. Six hundred million random density matrices were employed. 



V. DETERMINANTAL PRODUCT MOMENT FORMULAS 



A. Two-rebit case 



At a still later point in our investigation, we realized that we might make further progress- 
despite apparent limitations on the number of determinantal moments we could explicitly 
compute-by exploiting the evident pattern followed by our newly-found formulas (|3]) and 
Q-in particular, the structure in their denominators. This encouragingly proved to be the 
case, as we were able to additionally establish that 

<|Pl1p PT | 3 > 2 -^5=|<IH fc > 2 _ re ^ 5; (11) 

where 

A 3 = 8fc 9 +180fc 8 +1674A; 7 +8559A; 6 +29493A; 5 +84291A; 4 +136801fc 3 -401334A; 2 -2516616A;-3612816 

(12) 

and 

B 3 = 32768(fc + 3)(A; + 4)(A; + 5)(4A; + ll)(4A; + 13)(4A; + 15)(4A; + 17)(4A; + 19)(4fc + 21). (13) 
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FIG. 3: Difference obtained by subtracting the two-qubit HS probability distribution in Fig. [2]from 
the two-rebit probability distribution in Fig. [TJ Darker colors indicate more negative values. 

So, it then became rather evident that we can write for general non-negative integer n, 

A 



(\p\ k \p PT \ n ) 



2-rebit/HS ft 



2-rebit/HS 



(14) 



where both the numerator A n and the denominator B n are 3n-degree polynomials (thus, 
forming a "biproper rational function" [25]) in k (the leading coefficient of A n being 2 n ), 
and 



B n = 128 n (fc + 3) n 2k 



11 



(15) 



2n 



where the Pochhammer symbol (x) n = r jf^"^ = x(x + 1) . . . (x + n — 1) is employed. Further 
still, moving upward to the next level (n = 4), we determined that 

A, 



HV T IV 



2-rebit/HS D 



K> 



2-rebit/HS 



(16) 



where 



A 4 = lQk 12 + 576k 11 + 9112A; 10 + 84496A; 9 + 525681 A; 8 + 2389416A; 7 + 7805462A; 6 + 13904508fc 5 + 



(17) 



+6212189A; 4 + 166748972A; 3 + 1636873812A; 2 + 5496485760A; + 6610161600, 
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and -B4 is given by (15) with n = 4. The real part of one of the roots of A4 is 2.999905, 
suggesting to us some possible interesting asymptotic behavior of the roots of these numer- 
ators, n — > 00. In a related predecessor study [IS] [sec. II. B. 2], we had been able to discern 
the general structure that the denominators of certain "intermediate [rational] functions" 
used in computing the (univariate) moments of {p PT \ n ) 2 - re bit/ hs' ^ — 1> • • • j 9 followed. 

From our four new two-rebit determinantal moment results (|3]), Q, (11) and (16), we 
see that the constant terms in the 3n-degree numerator A n are —16,4860,-3612816 and 
6610161600 for n = 1,2,3,4. Since we had previously computed [19j[eqs, (33)-(41)] the 
moments of (|p PT | n ) 2-rebit/HS' ra = 1) • • • ) 9, we were also immediately able to determine the 
next five members of this sequence {—16, 4860, —3612816, 6610161600}. However, no general 
rule for this sequence, which would, interestingly, directly allow us to obtain a formula for 

(\p PT \ n ) 2 -rebit/HSi had ^ emer S ed fOT tllem - 

Certainly, it would be of interest to conduct analyses parallel to those reported above for 
metrics of quantum-information-theoretic interest other than the Hilbert-Schmidt, such as 
the Bures (minimal monotone) metric [U [TBI ES] • The computational challenges involved, 
however, might, at least in certain respects, be even more substantial. 

B. Use of Cholesky decomposition in rigorously finding formulas for general k 

After having posted the results above, along with additional ones, as a preprint [27], 
Charles Dunkl detailed a computational proposal that he had outlined to us somewhat ear- 
lier. The particularly attractive feature of this proposal was that it would-holding the expo- 
nent n of |p PT | fixed-be able to compute the adjustment factors for general k, rather than 
having to do so for sufficient numbers of individual members of the sequence k — 1, . . . , N, 
so that we could successfully apply the Mathematica command FindSequenceFunction, as 



had been our strategy heretofore. The proposal of Dunkl (Appendix D ) involved parame- 
terizing 4x4 density matrices in terms of their Cholesky decompositions. The parameters 
(ten in number for the two-rebit case and sixteen for the two-qubit case) would be viewed 
as points on the surface of a unit (due to the trace requirement) 10-sphere or 16-sphere. 
The squares of the points lie in a simplex. One can then employ the corresponding Dirich- 
let probability distributions over the simplices to determine the associated expected values 
(joint moments). (A further highly facilitating aspect here is that both |p| and the jacobian 
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for the transformation to Cholesky variables are simply monomials in the variables.) Using 
this approach, we were able to extend our single (n — l,a — 1) two-qubit result ^ to the 

n = 2 case, 

<HVT> 2 _^ 5 = ( 18 ) 
k{k{k{k{k{k + 15) + 67) + 45) + 220) + 4260) + 10944 



64(2fc + 9)(2fc + ll)(4fc + 17)(4Jfe + 19)(4fc + 21)(4Jfe + 23) 
Additionally, in the following array, 



2-qubit/HS ' 



V 



-16 


4860 


-3612816 


6610161600 


-23680812672000 


147885533254368000 


5 


2940 


-2516616 


5496485760 


-21644930613600 


144374531813568000 


9 


709 


-401334 


1636873812 


-7755993054000 


58524043784903280 


2 


368 


136801 


166748972 


-1199508017652 


11977854861441312 




203 


84291 


6212189 


-4378482660 


1052189083196640 




48 


29493 


13904508 


29246867605 


-30302414250528 




4 


8559 


7805462 


7876634465 


-6899036908859 






1674 


noon /i 1 c 

2389416 


2649513956 


3583820785224 






180 


525681 


883461210 


1632448582425 






8 


84496 


219916945 


477741210624 








9112 


40679505 


118164517947 








576 


5660714 


23817008856 








16 


575800 


3786901675 










40000 


469728096 










1680 


44685468 










32 


3143808 












153360 












4608 












64 



/ 



(19) 



we show (n — 1, ... ,6), column- by-column, the (3n+l) coefficients of the numerator polyno- 
mials in ascending order-the entries in the first row corresponding to the constant terms,. . . - 
in the two-rebit case. 

Additional results for the cases n = 7, . . . , 13 were found [23 eqs. (17)-(21)]. The 
leading (highest-order) coefficients in these thirteen sets of two-rebit results were found to 
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be expressible in descending order as 

CWi = 2 n ; C 3n = 3 x 2 n - l n{n + 2); CW-i = 2 n ~*n{n{n{9n + 32) + 24) - 45); (20) 

C 3n _ 2 = 2 n ~ 4 n (n (n (n (9n 2 + A2n + 52) - 119) - 52) - 60) . (21) 
From these four formulas, we are able to reconstruct (n = 1) all four entries in the first 



column of the table (19). Thus, it appears that, in general, C$ n -i is a polynomial in n of 
degree 2(z + 1). (For i = 3n — 1, we obtain the constant term, of strong interest. With the 
full knowledge of all the constant terms, and none of the other coefficients, we could obtain 
the univariate moments {\p PT \ n ) 2 -reUt/ hs'^ Further, we have found that 

C*3n-3 = (22) 



in— 7, 



5 

and 



- 1) (13577 7 + 855T7 6 + 189577 5 - 177177 4 - 3091?i 3 - 7731n 2 + 32394n) 



^ 3 n-4 = \2 n -\n - l)n (23) 



n(n(n(n(n(3n(3n(9n + 59) + 377) - 2887) - 2295) - 10535) + 112240) - 181492) + 436720. 
C. Two-qubit formulas 

The numerators of our four sets (n = 1,2,3,4) of two-qubit results (the first two having 
been obtained by "brute force" Mathematica computations, and the last two, using the 
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Cholesky-decomposition parameterization) are expressible, in similar fashion, as 



/ 



V 



-42 


10944 


-6929280 


9247219200 


-1 


4260 


-3684384 


6039653760 


6 


220 


-456948 


1342859616 


1 


45 


80168 


64072440 


— 


67 


27783 


-13235252 


— 


15 


5373 


1080858 


— 


1 


1458 


1160375 


— 


— 


282 


278478 






27 


50991 






1 


7542 








749 








42 








1 



\ 



(24) 



/ 



We observe that the leading coefficients C^ n+ \ of all four numerators are 1, so they are monic 
in character, while the next-to-leading coefficients fit the pattern C^ n = 3n(n + 3)/2. 



It is evident at this point, in striking analogy to the general two-rebit formula (14), that 
in the two-qubit scenario, 

A, 



(\p\V t \% 
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2-qubit/HS r> /2-qubit/HS 

where, again, both the numerator A n and the denominator B n are 3n-degree polynomials in 

9" 



k, and (cf. (15)) 
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VI. DETERMINANTAL PRODUCT MOMENT FORMULAS FOR 6x6 DENSITY 
MATRICES 

Of course, one may also consider issues analogous to those discussed above for bipartite 
quantum systems of higher dimensionality. To begin such a course of analysis, we have 
found for the generic real 6x6 ("rebit-retrit") density matrices (occupying a 20-dimensional 
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space) the result 

, k PT , Ak b + 40A; 4 + 95A; 3 - 220A; 2 - 1149Jfe - 1170 , fc . 

\\P\ IP \)rebit-retrit/HS ~ 57 Q( k + 4 )( 3fc + 11 )( 3fc + 13 )( 6fc + 2 3)(6k + 25) ^ ' rebit-retrit/ HS ' 

(27) 

Increasing the exponential parameter n from 1 to 2, we obtained that the rational function 
adjustment factor for (\p\ k \p PT \ 2 ) rebit _ retrit / HS is the ratio of 

16fc 9 +336fc 8 +2616fc 7 +8496fc 6 +12069fc 5 +101979A; 4 +903539A; 3 +3316809A; 2 +5620320A;+3715740 

(28) 

to another ninth-degree polynomial 

331776(A; + 5)(3A; + 11)(3A; + 13)(3A; + 14)(3A; + 16)(6A; + 23)(6A; + 25)(6A; + 29)(6A; + 31). (29) 

Additionally, for the generic complex 6x6 (qubit-qutrit) density matrices (occupying a 
35-dimensional space), we have obtained the result 

, k PT , _ k 5 + 15fc 4 + 37k 3 - 423A; 2 - 2558A; - 3840 , k , 

\\P\ \P \) qubit-qutrit/ HS ~ 72 {2k + 13)(3fc + 19)(3fc + 20)(6fc + 37) (6k + 41) ^' P ' / qubit-qutrit/ HS 

(30) 

It should be pointed out, however, that in contrast to the 4x4 density matrix case, the 
nonnegativity of the determinant of the corresponding partial transpose of a 6 x 6 density 
matrix does not guarantee separability, since possibly two eigenvalues of the partial transpose 
could be negative, indicative of entanglement, while still yielding a nonnegative determinant 

(cf. nz|). 



VII. MINIMALLY DEGENERATE TWO-REBIT DENSITY MATRICES 

For the eight-dimensional manifold composed of generic minimally degenerate two-rebit 
systems (corresponding to density matrices p with at least one eigenvalue zero), forming the 
boundary of the nine- dimensional manifold of generic two-rebit systems, we have computed 
the Hilbert-Schmidt moments of |p PT | n , n — 1, . . . , 10. (For such systems, |p PT | G [—jq, 432]-) 
These results are given in Appendix [Cj (Charles Dunkl was able to find rational functions of 
k for n = 1, 2, 3-but not yet further- which yielded these moments when k was set to zero.) 

We note that as a particular case of results of Szarek, Bengtsson and Zyczkowski |28j, 
the Hilbert-Schmidt probability that a generic two-rebit system is separable is twice the HS 
probability that a generic minimally degenerate two-rebit system is separable. 
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FIG. 4: Two sets of estimates of the Hilbert-Schmidt two-rebit separability probability. The 
upper (blue) decreasing curve is based on the first 3,310 (nonnegative) moments of (\p\ \p PT \) n 
and the lower (red) increasing curve on the first 3,310 (alternating in sign) moments of \p PT \ n . 
The true separability probability, thus, appears constrained to lie within the range [0.453104500, 
0.454543513]. 

VIII. ESTIMATION OF SEPARABILITY PROBABILITIES, USING CONJEC- 
TURED FORMULAS 

A. Two-rebit case (a = \) 



We now utilize the conjectured formulas (App. D 6 (-developed by Dunkl at an intermedi 



ate stage in our research effort-with the Dyson- index-type parameter a set to |, correspond- 
ing to the two-rebit case. In Fig. [4], we display the corresponding Hilbert-Schmidt separabil- 
ity probability estimates obtained by application of the Legendre-polynomial-based probabil- 
ity density reconstruction (Mathematica) procedure of Provost [29] ec L- (15)]— yielding least- 
squares approximating polynomials-to the sequence of the first 3,310 moments of (|p||p PT |) n 
(upper blue curve) and to the sequence of the first 3,310 moments of |p PT |"' (lower red curve). 
(All our computations here and below were conducted with 48-digit accuracy. A uniform 
"baseline density" was, in effect, assumed, while the use in this capacity of a beta distribu- 
tion, fitted to the first two moments, and Jacobi polynomials yielded highly erratic estimates 
when the corresponding Mathematica algorithm of Provost j29j pp. 750-752] was applied.) 

In Fig. [ij the last/highest pair of estimates is {0.453104500, 0.454543513}, so it certainly 
appears that the true (common) separability probability for the two variables must lie within 
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this interval. The convergence properties of the two sequences of estimates display paral- 



lel (increasing-decreasing) behavior in the two-qubit case. (In sec. D5, Dunkl develops a 
distinct/alternative probability distribution reconstruction approach of interest-which he 
applies to considerably fewer moments than the 3,310 we do-to the two-rebit separability 
probability estimation problem.) 

Our 2007 hypothesis ([21 sec. X.A]) that the Hilbert-Schmidt separability probability of 
generic two-rebit systems is ~ 0.470588 can, thus, be decisively rejected (Fig. [Ij), since 
it clearly lies outside the confining interval. We will here note that in the later 2010 study 
[7] [p. 7], a numerical estimate of 0.4528427, substantially different from ^, was reported, 
and it was additionally observed that in [21 sec. V.A.2] the best numerical estimate of the 
two-rebit separability probability obtained there had been 0.4538838. A possible exact value 
of || = 0.453125-which does lie within the confining interval in Fig. [IJ-was, in fact, suggested 
in p. 6]. Use of linear algebraic principles, did allow us in [7j to establish an upper bound 
on the generic two-rebit Hilbert-Schmidt separability probability of ~ 0.537619. 

We note, importantly, that the lower bound of the confining interval, 0.4531014500 is 
0.999955 times as large as jj| . 



B. Two-qubit case (a = 1) 

In Fig. [5] we similarly show-for the two-qubit case (a = l)-the estimates obtained by 
application of the probability distribution reconstruction procedure of Provost [221 eq. (15)] 
to sequences of 2,415 moments of (|p| |p PT |) n (upper blue curve) and |p PT | n (lower red curve). 
We, of course, note that the lower bound obtained of 0.2424235313 seems to nicely support 
our 2007 hypothesis ([21 sec. X.B]) that the Hilbert-Schmidt separability probability of 
generic two-qubit systems is ^ ~ 0.242424. (The ratio of this lower bound to that based 
on 2,414 moments is 1.000000006779, indicative of strong convergence. The analogous ratio 
for the upper estimate was 0.99999153401-somewhat less strong.) 

Zyczkowski, Horodecki, Sanpera and Lewenstein, in their foundational paper [H eq. (36)], 
provided a numerical estimate-0.632±0.002-of the generic two-qubit separability probability, 
using as a measure the product of the uniform distribution on the 3-simplex of eigenvalues 
and the Haar measure on the 15-dimensional 4x4 unitary matrices. (The 4x4 density 
matrices were, then, in a sense, over-parameterized. The authors were "surprised" that the 
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FIG. 5: Two sets of estimates of the Hilbert-Schmidt two-qubit separability probability based on 
2,415 moments. The upper (blue) decreasing curve is based on the (nearly all nonnegative) moments 
of (|/o||p PT |) n and the lower (red) increasing curve on the (alternating in sign) moments of \p PT \ n . 
The true separability probability, thus, appears to lie within the confining range [0.2424235313, 
0.2485026468].) 



probability exceeded 50%.) They also advanced [U, eq. (35)] certain analytical arguments 
that the probability was in the interval [0.302, 0.863]. While these studies are of great 
conceptual interest, they did not specifically employ as measures those denned by the volume 
elements of metrics of interest (such as the Hilbert-Schmidt, Bures,. . . ) over the quantum 
states. 



C. Reconstructed probability distributions 

In Figs. [6] and [7| we show (based on 200 moments, using now the procedure of Mnat- 
sakanov [30] ) , rather than that of Provost [29] , the reconstructed HS two-rebit and two-qubit 
probability distributions for both sets of moments, all distributions linearly transformed to 
the interval [0,1]. 



D. a as a free parameter 



As an exercise of interest, let us consider the Dyson-index-like parameter a in sec. |D 6 
with the values \ and 1 conjecturally corresponding to the two-rebit and two-qubit moments, 
respectively, as a free / 'continuous parameter (cf. [31]), and perform our standard separability 
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FIG. 6: Reconstructed-and linearly transformed to [0,1]-HS two-rebit probability distributions 
based on 200 moments of |/?||p PT | (blue, lower-peaked curve) and \p PT \ (red, higher-peaked curve) 
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FIG. 7: Reconstructed-and linearly transformed to [0,1]-HS two-qubit probability distributions 
based on 200 moments of |/?||p PT | (blue, lower-peaked curve) and \p PT \ (red, higher-peaked curve) 



probability calculations using the Provost algorithm [29j-taking the same ranges as before 
for the determinantal moment variables. Based on ninety-six moments, we obtain Fig. 18} 



E. a = 2 (quaternionic?) 



In Fig. M we show-for the a = 2 (presumptively quaternionic) case (Appendix D 6 )-the 



estimates obtained by application of the procedure of Provost j29j eq. (15)] to the sequences 
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FIG. 8: Separability probability estimates as a function of the parameter a (sec. D 6). The upper 



curve is based on ninety-six moments of \p\\p |, and the lower curve on ninety-six moments 



of \p |. Also included as horizontal lines are the two-rebit (a 
two-" quaterbit" (a = 2) and "classical" (a = 0) conjectures of 



2/' 



two-qubit (a = 1) and 
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0.242424, 
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323 



0.080495 and 1, respectively. 



of moments of (|p||p PT |) n (upper blue curve) and \p PT \ n (lower red curve). (We use the 
term "presumptively", precisely because we have performed no explicit calculations-as we 
certainly have done in the two-rebit (a = |) and two-qubit cases (a = l)-involving 4x4 
quaternionic density matrices. We are, thus, proceeding under the assumption that we 
can extrapolate the formula of Dunkl to the case a = 2. Dunkl, however, has noted that 
his formula does agree with that of AndaipiH Thm. 4], in the quaternionic case, for the 
[univariate] moments of \p\ (cf. [32]). Also, Dunkl has raised the issue of whether or not 
nonnegativity of the determinant of the partial transpose is equivalent to separability, as 
it is known to be in the two-rebit and two-qubit cases |17j.) The lower estimate based on 
2,325 moments is 0.080495355 (which is 1.000000000049 times the corresponding estimate 
based on 2,324 moments). This 2,325-moment estimate can be truly remarkably well-fitted 



26 
323 



0.0804953560. 



by the relatively simple fraction 

In the framework of [2] [sec. IX], the "scaling factor" used to obtain the ^ result would 
be SqS , where 19136 = 2 6 x 13 x 23 and 152809335 = 3 6 x5x7x 53 x 113. (In 
these calculations, we took the total HS quaternionic volume to be equal to the product of 
that volume given by Andai in [16] and the normalization factor of 2 13 indicated there-thus, 
giving us the HS volume in the Zyczkowski-Sommmers framework [3] that we have employed 
throughout.) For our two other conjectures, the associated scaling factors would be (a = |, 
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FIG. 9: Two sets of estimates of the (quaternionic?) Hilbert-Schmidt separability probability. 
The upper (blue) decreasing curve is based on Dunkl's conjectured formulas-using a = 2-for the 
expected values of and the lower (red) curve, similarly for \p PT \ n . 2,325 moments were 

employed. 
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cases, respectively. 



and (a = 1, two-qubit) . The associated HS separable volumes would, 
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, for the real, complex and quaternionic 



F. a = 4 (octonionic?) 



In Fig. 10 we show-for the a = 4 (octonionic? (cf. [9j QUI US])) case-the estimates 



obtained by application of the procedure of Provost j29j eq. (15)] to sequences of 2,125 
moments of (|p||p PT |) n (upper blue curve) and \p PT \ n (lower red curve). The fraction 



760 



2 3 -5-19 



69903 — 3 4. 863 ~ 0.0108722086 is 0.9999999981 times as large as the estimated separabil- 
ity probability. Convergence is comparatively very strong in this instance, and definitely 
seems to improve, in general, as the Dyson-index-like parameter a increases. 



G. a = (classical?) 



If we set a = in (D 6) for the (mixed- moments) case n = k, we obtain the simplification 

4096- n r(2n + l) 3 



{W\p pt \ n ) 



(2)271 (2) 



(31) 



in 
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FIG. 10: Two sets of estimates of the (octonionic?) Hilbert-Schmidt separability probability. 
The upper (blue) decreasing curve is based on Dunkl's conjectured formulas-using a = 4-for the 
expected values of (\p\ \p PT \) n and the lower (red) curve, similarly for \p PT \ n . The true value appears 
to be constrained to lie within [0.0108722086, 0.0264396063]. 2,125 moments were employed. 
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FIG. 11: Two sets of estimates of the (classical?) Hilbert-Schmidt separability probability, using 



a = 0. The upper (red) decreasing curve is based on Dunkl's conjectured formula (D6) for 
the expected values of \p PT \ n -&nd the lower (blue) curve, similarly for (| / o||/9 P:r |) n -given by the 



simplified formula (31). 1,650 moments were employed, with the last pair of estimates being 



{0.96238936,0.99445741}. 



In Fig. 11, we plot our standard pair of two estimates (although now the roles of upper 



and lower curves are reversed). It appears that there is convergence to 1, that is, a = 
corresponds, in some sense, to a classical scenario, in which no entanglement is present. In 
regard to setting a = 0, Dunkl commented that doing so "assigns measure zero to the off- 
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diagonal entries of the Cholesky factor. The determinant and PT-determinant are identical 
as far as the measure is concerned, and the probability distribution is the same as that of the 
product M2M4 on the simplex in 3-space (ti,t 2 ,t 3 > 0, t 4 = l — t 1 — t 2 —t 3 and t 4 > 0)." His 
"attempt to reconstruct the underlying probability distribution yields an inelegant integral 
of a hypergeometric series" . 



H. Other values of a 



We also have conducted Legendre-polynomial reconstruction analyses for a number of 
other values of a, which we summarize in the form (cf. Fig. [8j 
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{0.64744667,0.63955009} 
{0.34299437,0.32784144} 
{0.13756171,0.14950325} 
{0.029008076, 0.055230359} 



\ 



(32) 



The first two columns give the value of a and the number of moments employed, and the last, 
the confining interval for the associated separability probabilities, with the first value being 
based on the moments of \p PT \ and the second, on the moments of |p||p PT |. Convergence 
of the probability-distribution reconstruction algorithm, based on the moments of |p PT |, 
appears to greatly increase as a increases. (An extremely close fractional fit to the lower 



bound for a = 8 is 



81 



318407 



0.00025439139215.) 



I. Specialized lower-dimensional ("non-generic") cases 

In [21 sec. II. A], we considered classes of 4 x 4 real, complex and quaternionic density 
matrices, where-as usual-the diagonal entries were allowed to take values in the 3-simplex, 
but now five of the six pairs of off-diagonal entries were nullified, leaving only the (2,3) and 
(3,2)-pair as free. (The associated separability probabilities were found to be ||, | and ^.) 



Dunkl (App. D7) has now been able to prove formulas for the bivariate moments in these 



specialized scenarios. 
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FIG. 12: Probability distributions (33) and (34) over t = 2 s \p\ (t <G [0,1]). The Hilbert-Schmidt 



(red) curve dominates the Bures curve above t = 0.021702. 



IX. HILBERT-SCHMIDT AND BURES PROBABILITY DISTRIBUTIONS OVER 



In the course of this work, Charles Dunkl further communicated to us a result (following 
his joint work with K. Zyzckowski reported in [31], where "the machinery for producing 
densities from moments of Pochhammer type" was developed) giving the univariate proba- 
bility distribution over t G [0, 1] that reproduces the Hilbert-Schmidt moments of t = 2 8 |p|, 
where p is a generic two-rebit density matrix. (If we set n = in our general [bivariate] 
determinantal moment framework above, we obtain the [univariate] moments of \p\.) This 
probability distribution took the form (cf. [E)][eq. (4.3)]) 

63 / r ~ r „ n \ _ , / r ~ \ i r. 



(^Jl-y/t (-8t - 9Vi + 2) + 15tlog (Jl-Vi + 1^ - ~tlog(t) 



(33) 



(see Appendix D 2 below for further details). At the suggestion of the author, Dunkl was also 
able to derive, in similar fashion, the Bures metric [U|26] counterpart of this Hilbert-Schmidt 
result (33). It took the form (Appendix D 3[ ) 



4V V~t - t (2y/i + 13) + 3tt + l) + 2 (l2y/i + 3) sin" 1 (l - 2y/i) 



71 



(34) 



In Fig. 12 we display these two (Hilbert-Schmidt and Bures) probability distributions. 
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X. DISCUSSION 



A. Background 

A basic linear-algebraic criterion that a Hermitian matrix be nonnegative-definite, that 
is have all its eigenvalues nonnegative, is that all its principal minors be nonnegative. In 
[7J, we were able to implement this criterion, in part, making use of the 3x3 minors, 
establishing thereby that the Hilbert-Schmidt probability a generic two-rebit system is sep- 
arable is bounded above by ~ 0.537619. (The absolute separability probability of 
6928 29/2 2057r ~ 0.0348338 provided the best exact lower bound established in this specific 
setting [7j, it appeared. The set of absolutely separable two-qubit states are described in 
Figs. 1-5 in [35] (cf. [361438] ). No immediate application of the moment-based approach 
adopted in this study to the description of the absolutely separable states is apparent.) That 
study [7] was a continuation of a series of papers of ours (including [2j [3914^6] ) in which we 
examined the separability probability question-for the Hilbert-Schmidt as well as various 
monotone (such as the Bures) metrics-from a variety of mathematical perspectives, employ- 
ing a number of density-matrix parameterizations. A major motivation in undertaking the 
moment-related analyses reported above was to further sharpen our separability probabil- 
ity estimates, perhaps even being able to arrive at an estimate accurate to several decimal 
places, and possibly obtain thereby convincing evidence for a particular true value. 

Despite the considerable computational efforts expended in calculating high-order mo- 
ments, the goal of high accuracy nevertheless appeared remote-that is, until the appar- 
ent advances of Dunkl (Appendix D) that we have sought to subsequently exploit above. 
This somewhat pessimistic viewpoint had been based on a continuing series of attempts 
by us-using a wide variety of probability-density reconstruction methodologies-to isolate 
the two-rebit separability probability on the basis of the initially computed (limited num- 



ber of) thirteen moments. As an example (cf. sec. D5), use of the nonparametric proce- 
dure of Mnatsakanov [30], yielded HS generic two-rebit separability probability estimates of 
0.4582596, 0.42970496 and 0.40321291 based on the first eleventh, twelfth and thirteen mo- 
ments of |p PT | (sec. [A| , so, no convergence was apparent, at least, with these few moments. 
The corresponding estimates were 0.5414052, 0.3923661 and 0.4792091 based on eleventh, 
twelfth and thirteen moments of |p||p PT | (sec. LB]). Use of the first ten moments in a certain 
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maximum-entropy reconstruction methodology [UJ gave an estimate of 0.409858. Addition- 
ally, incorporation of the first twelve moments into an adaptive spline-based algorithm [4"S] 
gave 0.4502338. The semiparametric Legendre-polynomial-based reconstruction approach 
of Provost [29]-our chief computational procedure in the main body of this paper-gave esti- 
mates of 0.3856787 and 0.4846628 based on the first thirteen moments of |p PT | and |p||p PT |, 
respectively. 

We had, thus, before the general formula of Dunkl, encountered evident difficulties in 
ascertaining to high accuracy the values of separability probabilities. These difficulties, it 
seemed, perhaps manifested the iVP-hardness of the problem of distinguishing separable 
quantum states from entangled ones [4"9Tl5"T] . As possible evidence for such a contention, if 
one knew all the generic HS two-rebit moments of |p PT |, then presumably one could deter- 
mine the associated separability probability to arbitrarily high accuracy. But to know all 
these moments, it appeared that one would have to know an indefinitely large number of the 



functions C*3 n _j ((j20[)-( 23 )), from which the needed constant terms could be extracted. In 
the apparent absence of a generating rule for these increasingly high-order functions (but see 
Appendix [D]), an indefinitely large amount of computation appeared to be required. ("Al- 
though [quantum entanglement] is usually fragile to the environment, it is robust against 
conceptual and mathematical tools, the task of which is to decipher its rich structure" [52l p. 
865].) "In [19j sec. II. B], an earlier study of ours of the moments for two-rebit systems, we 
encountered a somewhat analogous rather intractable state-of-affairs, employing the Bloore 
(correlation-coefficient) parameterization of density matrices (and not the Cholesky decom- 
position parameterization, as in this study). There, a general formula for the denominators 
of certain important "intermediate functions" could be discerned, but only explicit results 
obtained for an initial set (m — 2,4,6, ... , 16) of the corresponding numerators. So, higher- 
order moments-and, thus, high accuracy-appeared out of reach there (but certainly in light 
of the apparent progress-but not yet rigorously established-of Dunkl, the matters there 
might also be readdressed) . 

B. Results 

In this paper, we have advanced four specific conjectures (a = 0, \, 1,2) (Fig. J8J). The 
reader might have been somewhat skeptical of our strong predisposition to conjecture ra- 
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tional values for the various separability probabilities under consideration. A basis for this 
inclination had been established in [2J, where a pattern of rational separability probabilities 
appeared through the application of exact methods to lower-dimensional non-generic (but 



more easily computed) quantum scenarios (sec. Villi). 

In regard to the conjecture (2J, sec. IX. B] that the Hilbert-Schmidt separability probability 
of generic (15-dimensional) two-qubit systems is K. Zyczkowski informally wrote: "It 
would be amazing if such a simple number occurs to be true! I wonder then if it is likely that 
this result may be derived analytically (by a clever integration), or perhaps even 'guessed' 
from some symmetry arguments [which are still missing]". From the author's viewpoint, 
perhaps one of the chief hurdles here is simply the exceptionally high-dimensionality and 
quartic (separability) constraints that need to be addressed in any integration ("clever" or 
otherwise). Possibly with the advent of more powerful symbolic (quantum?) computational 
systems, this obstacle might be directly overcome. Also, in terms of symmetry principles, the 
(Keplerian) concept of " Stella octangula" [Ml El] has proved useful in studying separability, 
and might conceivably do so (in some higher-dimensional realization) in the future. Certain 
interesting aspects of convexity were applied in (28] to obtain theorems pertaining to Hilbert- 
Schmidt separability probabilities. 

The general formulas of Dunkl remain formally unproven. However, our confidence 
in their validity is certainly enhanced by the reasonableness and non-anomalous behav- 
ior (Figs. [i| [5j |6j [7J [9j 10) of our various (separability) probability estimation procedures, 
for various values of a, which rely upon them. If the formulas did not, in fact, yield genuine 
moments of probability distributions, we would certainly expect that to be manifested, in 
some overt manner (negative probabilities, probabilities greater than unity, non- convergent 
behavior,. . . ) in our reconstruction efforts. 

It is interesting to note that of our three basic (two-rebit, -qubit, -"quaterbit" [TT] ) 
separability probability conjectures-||, ~-th.e two-qubit is the simplest, in the sense of 
having the smallest denominator (and numerator). The two-qubit systems exist conceptually 
in the framework of (standard/conventional/phenomenological) complex quantum mechanics 
P21E21 sec. 2]. 

A further observation is that although in random matrix theory, a (Dyson-index) param- 
eter P = 1 (the dimension of the corresponding division algebra [TU]) is typically assigned 
to the real systems, in the (Cholesky decomposition-based) analysis of Dunkl (App. [Pj), 
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the use, instead, of a = | appears to be natural-since one-halves repeatedly arise in the 
integration over the real sphere in TZ 10 . 

Knowledge of all the moments of |p PT | and |p||p PT | theoretically determines the complete 
probability distributions of these two variables (since the ranges of these two variables are 
bounded). In some sense, this constitutes more information than it might seem one should 
require to determine the single (separability) probability of primary, motivational interest 
[I]. So, if at some point in time, the separability probability questions can be resolved by 
some more direct methods, than it may appear that the analytical moment-based approach 
pursued here was more than was, in fact, truly required for the task at hand. Nevertheless, 
in the interim, this approach has clearly greatly advanced our knowledge of the ranges within 
which the separability probabilities must lie-even if not helping to pinpoint their conjectured 
exact (simple rational) values. 

C. Bures analyses 

In a naive exercise, we investigated whether or not the bivariate moment formulas pre- 
sented here might further hold-at least up to proportionality-if one were to simply replace 
the expectation with respect to the Hilbert-Schmidt metric in them by expectation with 
respect to the Bures (minimal monotone) metric [U [TBI ESI ESI E2] ■ However, such a possi- 
ble relationship appeared to be quite emphatically ruled out, at least with the one specific 
example, formula ^ above, we numerically studied in these regards. 

In [391 ec L- (16)] we had-based on extensive quasi-Monte Carlo numerical integrations- 
advanced the hypothesis that the two-qubit Bures separability probability took the form 
(with the "silver mean", o^g = V% — 1) 

PbLs = ~ 0.07333893767 (35) 

(which we do note is obviously irrational-in contrast to our Hilbert-Schmidt conjectures). 
We have recently begun to reexamine the results of that 2005 study, particularly in light of 
the later (2009) development, making use of Ginibre ensembles, of a "simple and efficient al- 
gorithm to generate at random, density matrices distributed according to the Bures measure" 
[55] (cf. [56l eq. (22)]). In an ongoing calculation, employing extended-precision indepen- 
dent normal random variables, we have obtained (using the normal approximation to the 
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binomial distribution)-based on 281,350,000 realizations (20,627,508 being separable, giving 
a probability of 0.0733162)-a 95% confidence interval {0.07328572,0.07334664}. We note 



that this interval does contain the conjectured value (35) for the true Bures two-qubit sepa- 



rability probability. (Consistently with these analyses, if we introduce our Hilbert- Schmidt 
two-qubit separability-probability conjecture of J| into the inequality of Ye [571 mid. p. 
7], we obtain 0.00373882 as a lower bound on the Bures two-qubit separability probability. 
Application of the very next inequality of Ye appears to yield 599089., obviously greater 
than 1, as an upper bound on this probability.) 



Appendix A: Two-rebit Hilbert-Schmidt moments (\p 



PT\n\ 



, 2-rebit/HS'> 



n 



13 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 



j_ 

858 

27 
2489344 

8363 



66216550400 

21859 
10443295948800 

23071 

539633583390720 

3317321 
3253917653076541440 

419856257 

15366774022001834065920 

16945249 

21117403549591928832000 

6102620963 

240565904621616585139814400 

87816716413 

103068223454742370906999357440 

7685831825319 

255310031843279606667374504181760 

23559692226221 

21217623285399369347467109090721792 

31283325154283 

736092406055063912488279599166259200 



-0.0011655 
0.0000108462 
-1.26298 x 10~ 7 

2.09311 x 10~ 9 
-4.27531 x 10" 11 

1.01949 x 10~ 12 
-2.73223 x 10~ 14 

8.02431 x 10~ 16 
-2.53678 x 10~ 17 

8.52025 x 10~ 19 
-3.01039 x 10~ 20 

1.11038 x 10~ 21 
-4.24992 x 10~ 23 



(Al) 
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Appendix B: Two-rebit Hilbert-Schmidt moments {{\p\\p PT \) n ) 2 _ re ^/^5? n = 1, . . . , 13 



2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 



5696343244800 

1 

677899511057612800 

1 

45973294808920227840000 

1 

11662680803407302839532257280 

3929 

4158654163938276392103553381781471232 

1 

158158366213274948625327048295175946240 

71527 

1091771390479438557169317171313498708778365747200 

4847 

8524774835462825812953111833131999123882778862551040 

2637 

441859421690475898778224458156196857558486829112995348480 

1 

16833241044745336849504728327369136893812113975649318731776 

66838003 

103562821755098721107694750210986399334006111231977898090691403395891200 

55601 

7991978474394124344137676945763648296251134760058623476703807122125619200 



0. 

1.22886 x 1(T 12 
1.47514 x 1(T 18 
2.17518 x 1(T 23 
8.57436 x 1(T 29 
9.44777 x 1(T 34 
6.32278 x 10 
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6.55146 x 10~ 44 
5.68578 x 10- 49 
5.96796 x 10~ 54 
5.94063 x 10~ 59 
6.45386 x 10~ 64 
6.9571 x 10~ 69 



(Bl) 



Appendix C: Moments of \p 



PT\n 



n 



1, . . . , 10, for minimally degenerate pairs of 



rebits 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



2376 
7 

380160 
9 

34777600 

443 
89942261760 

461 

4032782401536 

5455 

1785064543223808 

631 

6948198442598400 

474017 

161763811601154048000 

4003573 

39645007353595350220800 

3397 

924892257224239349760 



-0.00210438 \ 
0.0000184133 
-2.58787 x 10~ 7 

4.92538 x 10" 9 
-1.14313 x 10~ 10 

3.05591 x 10~ 12 
-9.08149 x 10~ 14 

2.9303 x 10~ 15 
-1.00986 x 10~ 16 

3.67286 x 10- 18 j 



(CI) 
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Appendix D: Two-rebit and two-qubit moments 



Charles F. Dunkl[H][HI] 

Let Q denote the set of 4-by-4 (symmetric) real positive definite matrices, and let fli 
denote the matrices of trace one in Q. Recall (X) denotes the expectation of the random 
variable X, with the associated probability density being implicit from the text. Furthermore 
\p\ denotes det p. 

1. Construction of density functions 

We describe the tools used to determine densities whose moment sequence is given in 
Pochhammer form. Here we restrict to densities supported on [0, 1]. Let / (x) be defined on 

< x < 1, such that / (x) > 0, / is continuous on < x < 1 and f (x) dx = 1. There is 
an associated random variable X, with Pr {a < X < b} = f (x) dx. The moment sequence 
is (X n ) = Jg 1 x n f (x) dx, n = 0, 1, 2 . . .. Observe that the moment sequence uniquely defines 
the density because the support is a bounded interval. 

First we consider a beta-type distribution: let a, /3 > 0, and 

/ (x) = J77~~o\ xa l (1 " ^ , < x < 1, (Dl) 

[ X n f(x)dx= ( " } " ,71 = 0,1,2,.... 
Jo \ a + P)n 

(Recall B (a, (3) = ^v^gy •) This uses the identity T (a + n) /T (a) = (a) n := 
nr=i ( a + * — 1)> ^ ne Pochhammer symbol. 

Lemma D.l Suppose X±,X2 are independent random variables on [0,1] with densities fa, 

1 = 1, 2. Then the density for XiX 2 is 

If the moments of Xi,X 2 are pn = {X^) = Jq x n fi(x) dx then (X™X%) = pn pn\n = 
0, 1, 2 . . .. , that is, 

[ x n f(x)dx = pUp%\n = 0,1,2,.... 
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The Lemma was stated and used in [3U p. 123521-20]. Also we use the duplication 
formulae for Pochhammer symbols: 



(«) 2 „ 




a- 


{2n)\ 


= (l) 2 n = 


2 2n n! 


(2n + l)! 


= (2) 2n = 


2 2n n! 



2,' ' 

3 - 



2. Density of the determinant under the Hilbert-Schmidt metric 

The 10-dimensional cone fl is equipped with the measure rii<j<j<4 dpij (where p = 
(Pij)tj=i i s the generic matrix). The probability distribution on Q\ is the (9-dimensional) 
restriction of this measure. 

The following lemma applies to N-by-N positive-definite matrices for any N = 2,3, . . .. 
Each such matrix p has a Cholesky decomposition: 

P = c l c, 

where C is upper triangular with entries cy, q,- = for i > j and Cu > for all i. The 



entries of p are pij = Y^k=i c kiCkj = °kiCkj- Consider the Jacobian matrix |^ where 

the dependent variables are p^, i < j. 

Lemma D.2 Suppose p = C l C then 



det 



dp 
dc 



N 



2" IT 



JV+1- 



i=i 



Proof We use the simple fact: suppose yi = f\ (xi, x 2 , ■ ■ ■ , 1 < i < iV then the matrix 
(j^r j is lower-triangular (0 for j > i) and det (j^- j = YliLi Now order the (inde- 
pendent) variables: c n , c 12 , ... , c 1N , c 22 , • • • , c 2 n, c 33 , . . . cjv-i,jv-i, c N - ljN , c NN . For % < j, 
Pij - J2l=\ c kiCkj + cuCijcind thus 

N N n N 



det 



dp 
dc 



iV iV r\ iV 

nn| 2 =n>.Y- +i )- 

1 = 1 1=1 J 1=1 



Now set N = 4. The pre-image S of Q± (for the map C h-> C*C) is a modified octant of 



the unit sphere in IR 10 , because Tr (C*C) = 



i<«<i<4 u «i 



c 2 „-. Recall the condition c 7) - > 0, but 
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the other entries can have arbitrary signs. The surface measure dm (C) on S is essentially 
a Dirichlet measure: consider a monomial on S, that is, 



/(c) : = n 

i<*<i<4 



C ■ ■ 



then 

1. if rijj is odd for some i < j then J* s / (C) dm (C) = 0, 

2. if riij is even for each i < j then 

/ /(0)(MC)= F7 —^ — ^ n rg + ^ 

3. if Tijj is even for each % < j, and iV := Xa<j<j<4 n *i then 

/ /(C7)dm(C7) = 7 i- (l) ■ 

In our usage either case 1 or case 3 applies. Combining the Jacobian and the fact |p| 
c ii c 22 c 33 c 44 we obtain (for normalized measure, that is 7 r 5 detgrfm(C) = 1), k = 0,1,2,. 

1. if riij is odd for some i < j then J 5 \p\ k f (C) det §^dm (C) = 0, 

2. if riij is even for each i < j, and iV := Xa<i<j<4 n *i then 

7 / \ P \ k f(C) det %m(C) = —± 



Ak+N/2 



X 



(9) ( 2 )fc+n 22 /2 ( 9 ) ( 1 )fe+n 44 /2 II (9) 

\^/k+rau/2 V^/ fc+W2 i V ^ / 



fe+nn/2 \ / k+nxi/2 l<j<j<4 ^ ' "-ij/ 2 

The special case / (C) = 1 provides the moments of the random variable (p) ; indeed 



f I \k A + ^ ^ (l)fc( 2 )fc (DjJ 1 )* 



We know the range of \p\ is[0, ^g] (the maximum is achieved at p = to use the previous 
results consider X = 2 8 \p\. Then 

(l)„(2)n(l) n (l)n_ o4 n2- 2 "(4) 2n 2- 2 ™(2) 2n 



(X n ) = 2 



8n 



2 4 ™ ( 5 )2n ( 1 2 1 ) 2 n ( 5 k ("^n 

( 5 k (t) 2 « 
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Thus X is (equidistributed as) the product of two independent random variables Xi,X 2 
with 



(X{ 



n\ _ ( V2n 



(5) 2n 4 + 2n 2 + n' 



(2) 2n 



(-) ' 

V 2 / 2n 

Clearly X x has the density f\ (t) = 2t, < t < 1. The density of X 2 is 

1 / \ 5/2 

because 

The density / (x) of X is given by 



(is 
s 

5/2 C?S 



25 (2, 1) 



2t ^ 



^ (2, 1) 

G3 ' r /' tT 3 (l-«) 



2 -'V^ 



u~ 3 (1 - uf' 2 du 
5/2 du. 



The integral is evaluated as follows: set u — 1 — s 2 , rf-u = —2sds, 

/ (x) = 63a: / -^ds 

Jo (1 - s 2 ) 

63x J -s (15 - 25s 2 + 8s 4 ) 15 (1 + s) 



H In 



s=0 



8 | (1-s 2 ) 2 2 1-s 2 

= y |(1 - v^) 1/2 (2 - 9-v/x - 8x) + 15x In ^1 + ^ 1 - ^ - y^lnxj . 

Also / (x) — O ^(1 — a;) 7 ^ 2 j near x — 1. 
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3. Density of the determinant under the Bures metric 



Using the Bures metric one obtains 

(2)^1 (-On (2)20 



-8n 



") 



(1) 



2d 



2 8n (!) n (2) n (4) 2n (n+l)(2n+l) (4) 2n ' 
for n = 0, 1, 2, . . .. As above we consider the random variable X = 2 8 
The density / (x) of X, for < x < 1, satisfies 



1 



(!) 



2d 



(X) ^ (n + 1) (2n + 1) (4) 2n 



n = 0,1,2, 



We express X as the product of two random variables. 
Let 

A(t) = r 1/2 -i,o<t<i, 

then 



t n fi(t)dt 



1 



Next observe (from equation Dl): 



r(4) 



(n + l)(2n + l)' n "°' 1 ' 2 '---- 



r (!) r (§) io " " ^ ~ (4)„ 

so set s = t 1//2 (and note p^ws^ = ^ s = \t l ^ 2 dt) to obtain 



Vr 1/4 fl -t 1/2 ) 3/2 dt = 

t ./o 1 j (4) 2n 



2d 



n = 0,1,2.... 



Let 



/ 2 (t) = -r 1 / 4 (i-t 1 / 2 ) 3/2 ,o<t< 1. 

7T 



By Lemma D.l the desired density function is 

- 1 ■ ■■■ dt 



/(*) 



7Y 



A(f)AW 



TT\/X 



{t 1 ' 2 - x 1 ' 2 )^ (I -t^f' 2 dt. 



Substitute t = s 2 then 



/(*) 
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an elementary integral; indeed 



/ (x) = |3tt (40 + 1) - 4 (13 + 20) \Jy/i-x - 2 (120; + 3) arcsin (20 - l) j . 
As with the Hilbert Schmidt metric, / (x) — O ^(1 — a;) 7 ^ 2 j near x = 1. 



4. The joint moments of |p| and |p PT | 

The partial transpose p PT of p is obtained by interchanging the values of pu and p 2 3 (and 
p 4 i and P32). In this section we introduce a conjecture for 

p PT | n |p| fe ),A;,n = 0,l,2,3,..., 

using the density on Qi coming from the Hilbert-Schmidt metric. 
For the upper triangular matrix C and p = C l C we find 

\p PT \ 

= C 11 C 22 C 33 C 44 + 2CnC2 2 (CllCi4 - Ci 2 C 13 - C22C23) 
X (-C11C23C44 - C23C34C11 + C11C33C34C24 + C 22 C3 3 Ci4 

- C22C13C33C34 - C12C33C24 + C 23 C33C3 4 Ci2) 

— (cuCu — C12C13 — C22C23) 2 

X (4c 22 C23CiiCi4 + C n C 44 + C n C 34 + C n C 24 - 2cnCi3C2 2 C 2 4 - 2C11C12C33C34 

22 2 2 2 2 2 2 \ 

- 2C11C12C23C24 + C 22 Ci 3 - 2C12C13C22C23 + C22C33 + C 12 C 33 + C 12 C 23 ) 

— (cnCi4 — C 12 C 13 — C 2 2C23) 4 . 

Of course I p I — ^11 ^22 ^33^44 * 

We introduce some utility functions (throughout n,k — 
0, 1,2,...). For a rational function F (k) = of define the degree to be deg (p) — deg (g). 

F (k) = (\p\ k ) = ^ )k ® k ® k ® k 



(10)« 

F 1 (n,Jb) = (|^r|p|*)/(| A ,|* 



F 2 (n,A;) = (|p| fc (|p^|-|p|)")/(H fc ), 

(A; + 1)„ (k + 1) (k + 2)„ (k+ f) 
i2 (n, fc) = F (n + fc) /F (fc) = ^ >nK (4 ^\ 0) ^ ^ 
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(Note F 2 (0, fc) — 1 = R (0, k)). The goal is to find (and prove) a closed form for Fi (n, fc), 
that is a general formula. Direct computation for n — 1, 2, 3 shows that F x (n, fc) is rational 
in k of degree 0; and at first glance, does not have an obvious formula (for the numerator). 
Some experimentation leads to the observation that F\ (n, k) — R (n, k) is of degree —2 
(verified only for small n). This motivates the investigation of the decomposition 



I PT\ n f H \ I \ n ~3 l\ PT\ I |V 

3=0 u/ 

.7=0 V/ 



F 2 (j, k + n- j)R(n- j, k) . 



For n = 1 we compute 



16 {4k + 13) (A; + 3)' 

this is an encouraging result, and it implies F 2 (1, fc) = — i e ^k+i3)(k+3) ' °^ degree —2. From 
the known value of F\ (2, fc) and the equation 

F 1 (2, k) = R (2, fc) + 2F 2 (1, fc + 1) R(l, k) + F 2 (2, k) 

we find 

F J( 2,*)- (* + 12)(2* + 7) 



256 (A; + 3) (k + 4) (4A; + 11) (4k + 13) (4A: + 17) ' 
This is of degree —3, rather than the hoped-for —4, and the factor (k + 12) is not of 
the "good" type, a divisor of (k + 1) 4 . So we try to modify F 2 (2, k) by adding a bit of 
F 2 (l,fc + 1)^(1, k); in fact 

F 2 (2, k) + j-^—F 2 (1, k + l)R (1, A;) 



fc + 1 v ' ' v ' 7 128 (fc + 3) (fc + 4) (4A; + 11) (4A; + 17)' 

of degree —4. We now have a "good" expansion of F x (2, k), namely 

2k (2 \ 

R (2, k) + — F 2 (l,k + l)R (1, fc) + ( F 2 (2, fc) + — F 2 (1, fc + 1) i? (1, fc) J . 

The terms are of degree 0,-2,-4 and each is an expression in linear factors. Next we 
consider F 2 (3, fc). This turns out to be of degree —5 (rather than —6). Some effort leads to 
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the satisfactory result: 

F 2 (3, k) + ^F 2 (2, k + 1) R (1, k) + (fc + 1 |( fc + 2) ^2 (1, fe + 2) i2 (2, fc) 

45 

~ ~2048 (fc + 3) (k + 4) (fc + 5) (4fc + 11) (4fc + 13) (4fc + 21) ' 

F 2 (2, t + 1) R (1, t) + (fc ^ ) - 1 | 2) f 2 (1, * + 2) R (2, t) 

9 (fc - 1) (k + 2) (2k + 3) 
~ 4096 (k + 3)(k + 4) (fc + 5) (4fc + 11) (4fc + 13) (4fc + 15) (4fc + 21) ' 

At this point there are enough examples to try to fit a formula to these expansions. Indeed, 
for < j < n let 

Cj (n, k) 



2°»(* + 3)„(2* + £) ai 



then 

n 

F 1 (n,k) = J2^(^k), (D2) 

is the conjectured formula. The degree of Cj (n, k) is —2j. We use the descending Pochham- 
mer symbol (o)( n ) = FKLi (° + 1 — «) = (— l) n (— o) n - If the conjecture is valid then for 
1 < j < n 

Cj (n, k) - 



j\ (n + k — 1) 

J " 1 

x S 7 ( n + fc ~ 1 )(j-i-<) F 2 (j -i,n + k + i-j)R(n + i- j, k) . 



'(2j-l) 



2! 
i=0 



This generalizes the examples found above. For generic k there is the expression 

(k + l) n (k + |) (Jfe + 2)„ 

Fi n, fc) = / w \ 11X ; " (D3) 

U j 26«(A; + 3) n (2A;+f) 2n 1 ' 



, -n,l,2, k, 2k 2n I , 

x -,/ : ; I - , , : 1 



This sum is a terminating balanced hypergeometric series, ("balanced" means the sum of 
the numerator parameters + 1 equals the sum of the denominator parameters.) However 
the 5F4-SUH1 is symmetric in (n, k) and the summation range is < j < min (n, k). When 
< k < n this omits the terms in the first formula for the range < n — j < ZL =^. For this 
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case the best way is to use equation (D2) (or else use (D3) with generic k to compute the 
rational function, then substitute the desired integer value for k). 
The special case k = is: 

Fr „ m _ 2(2n + l)! (2n)!(-2n-j) w /_«=>, -«=i >2l § _ 

11 j 2-(n + 2)(^) 2n + 2-(3)„m 2n 4 3 U-n,l-n,l + n' 

Another interesting special case is k = n: 

(\p PT \ n \p\ n )=F 1 (n,n)F (n) 
( 2 ™)<(§L -n,l,§,-4n-| x 



2l2n (T) 4 n(^ + 1 ) V-2n-l,-2n-|,|-n' 
The conjecture for iq (n, fc) has been checked by computer-aided symbolic algebra up to 
n = 13. 



5. Gaussian quadrature 

The method of Gaussian quadrature based on orthogonal polynomials can be applied 
to the density problem (see jiHl Thm. 3.4.2, p.48]). Suppose \i is a probability mea- 
sure supported on a bounded interval [a, b] and the moments are fij := J^x^dp (x). 
The orthogonal polynomials {P n (x) : n — 0, 1, 2, . . .} for /i (where P ra is of degree n and 
fVP n (re) d/x (x) = for < j < n) are determined by the moment sequence. Solve the 
linear system 

n-l 

a^i+j = -Atj+n, < j < ra - 1 

8=0 

to obtain the coefficients {a^} for the monic orthogonal polynomial 

n-l 

P„ (x) = x n + ajg*. 

2=0 

Then P n has n distinct zeros Ai < A2 < • . . < A n , contained in (a, b) .The structural constant 
h n = j^P n (x) 2 dp (x) = Y^=o Q-ifJ'i+n + A i 2n- The Gaussian quadrature rule with n nodes is 

n 

Gn (P) = y Wn,iP(Ai) , 

8=1 
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Then Q n (p) = j pdfi for all polynomials p of degree < 2n — 1. The sequence of discrete 
measures {Q n : n = 2, 3, . . .} converges weak-* (in the dual space of C [a, b}) to the mea- 



mulative distribution function of p. (Consider this as a sort of mid-point integration rule.) 
The orthonormal polynomials satisfy the three-term recurrence 



The most common approach to the computations is to find the coefficients directly 
from the moments. This is known to be a numerically ill-conditioned problem, so a relatively 
large number of significant digits must be used in the calculation. The algorithm of [591 
p. 476] with 30-digit floating-point arithmetic was used here. The computed values were 
checked for accuracy by evaluating the errors 



For the moments of 16 p and n = 20 we obtain 

A -.9501 -.9081 -.8587 -.8024 -.7402 
w .2714 • 10~ 10 .1397 • HT 8 .2416 • HT 7 .2337 ■ 10~ 6 .1553 • 10~ 5 

A -.6734 -.6032 -.5309 -.4581 -.3860 
w .7908 • lO" 5 .3293 -10" 4 .1171 - 10" 3 .3669 • 10" 3 .1034 • 10~ 2 

A -.3160 -.2495 -.1877 -.1317 -.08248 
w .2671 • 10~ 2 .6408 • 10~ 2 .1446 • 10" 1 .3111 - 10 1 .6499 • 10" 1 

A -.04104 -.008293 .01040 .02973 .04698 
w .1372 .3467 .3440 .4894 ■ 10" 1 .1994 ■ 10" 2 
One observes that a majority of the zeros are in [—1, 0] and most of the mass is contained 
in [—0.05,0.05]. With n = 30 we find 4 zeros in (0, ij; linear interpolation of the c.d.f. 
gives Pr { |p PT | > 0} ~ 0.42924. 

The distribution of 2 16 |p| |p PT | is somewhat more spread out over the interval. For n = 30 
we find 17 zeros in (0, 1), and linear interpolation yields Pr {|p PT | > 0} ~ 0.46129. 



sure p. The piecewise linear graph formed by consecutively joining [a, 0], (Ai + A2) , wi\ , 
W (^2 + A3) , Wi + w 2 ] , • • • \ (Aj + A i+ i) , Yjj=i w j > • • • > Ify 1] i s an approximation to the cu- 



xp n (x) = a n p n+1 (x) + (3 n p n (x) + a, n ^xVn-i {x) ,P-i = 0,pi = 1 



n 




1=1 
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6. Conjectures for the complex case 



Here we consider the question of moments of \q pt \ when p is a 4-by-4 Hermitian positive- 
definite matrix of trace one. The conjectured formulae have been verified for n = 1,2,3,4 
(see the previous sections). The conjecture was arrived at by inspecting the real case and 
using an analogous approach to the computed examples. It is interesting that the real and 
complex conjectured formulae can be combined into one formula with a parameter a. Set 
a = \ for the real case, a = 1 for the complex case. One could speculate whether a = 2 is 
related to a quaternionic or symplectic version [TBI E21 E3]- The general formulae are 

,fc\ _ k\(a + l) k (2a + l) k 



IP 



2« (3a + |), (6a + f) 

/(!/) 



2/v 



PT \ n i fe \ 1 
P P 



2 6 ™ (A; + 3a+ |) n (2£; + 6a + §), 

x (-2/e - 2n - 1 - 5a) J (/c + 1 + a) n ^. (/e + 2 + a)^. . 
For generic k this formula can be written as 

KTipi fc )/(H fc ) 

_ (* + l) n (k + 1 + a) n (k + 1 + 2a)„ 



2^ (k + 3a + l) n (2k + 6a + l) 2n 

—n, —k, a, a + |, — 2/c — 2n — 1 — 5a 
-k - n - a, -k - n - 2a, -*±* -*±n=l ' 



The special case n = is 



<ipr|p PT D 



(2n)! (1 + a) 2n (1 + 2a) 2n ( -n, a, a + \, -An - 1 - ha 



2l2n (3« + f ) 2n (6a + |) 4n ^ V -2n - a, -In - 2a, \ - n ] 1 
For k = we have 

n!(a + lW2a + l) n 



<ip pT r> 



26n (3« + i) B (6a+|) 



2d 



, (-2n-l-5a) w (a) B (a+|) n ^ / -==1, -n, a + 1, 2a + 1 \ 

2 4n (3a+ |) n (6a+ |) 2n 5 4 \1 - n, n + 2 + 5a, 1 - n - a, \ - n - a' /' 

because of the denominator parameter 1 — n it is necessary to replace the 5.F4-SUH1 by 1 to 

obtain the correct value when n — 1. 
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7. Lower-dimensional ("non-generic") case study 



In the Cholesky method, set five of the off-diagonal entries to zero; the positive matrix p 

o' 



is 



x\ 
x\ 



(where Xi > 0, 1 < i < 4, and x 5 comes from M. 13 , equipped with 



x 2 x 5 

x 2 x* 5 x\ + £5X5 

x\ 

an algebra structure including conjugation and a norm, e.g. (3 = 2, C ) 

Then \p\ = x\x\x\x\, and \p PT \ = x\ [x\ + x^x^) {x\x\ — x^x^x^) . Consider p as an 

element of M 4+/3 ; then the Jacobian for the map C ^ C*C is J = 16 Xl x% +1 x 3 x 4 . Write 

x$xl — \X5\ . With the usual Dirichlet integral techniques, integrating over the unit sphere 

Ylt=i + \x 5 \ 2 = 1 we get the normalized integral 



J x 2 ^xl m2 xl mi xl mA \x b \ 2ms (\p\) k J (x) dp 

(k + l) mi (k + 1) TO3 (k + l) m4 (k + l + l) m2 (l) m5 



(4 + (3 + 4k) 



S(k) 



\m\ 



where \m\ = Y^=i m ii an d 



*(*) := J (\p\) k J(x)dp 

{\p PT \) n {\p\fj{x)dp 

(k + l)l(k + l + l)l 
(4 + + 4k), n 

-n,k + l + £ + n,k + l + £ + n, § 
—A; — n, —A; — n, /c + 1 + | 

|2\ n 



^ 3 (l + f); 
(4 + 



Then 



Proof: Expanding 




f\ PT\\ n 2n I 2 1 I |2\ n / 2 2 2 | |2\ 

(IP I) = X 2 {X 3 + \X 5 \) (X 1 X 4 - X 2 \X 5 \ ) 



ij=0 



n\ In 



i 2n-2j' 2n+2j' 2n-2i 2n-2j 



X/1 



F5 



2i+2j 
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Now integrate with the above formula (value divided by 5 (k)) to obtain 
1 

x( * + 1) -(*L 



W:SG") ( - i),(4+i ^(* +i+ 9. 



by the Chu-Vandermonde sum, the second line equals (k + 1 + f + j) . Use the substitu- 
tions 



2 J A" (* + i + f), 



in the j-sum to produce the stated formula. 
Example: 



/ (|p PT |)(H) fe JW^ 



(4 + /3 + 4fc), 



x i (2fc + 2 + /5) | (k + l) 2 (2fc + 2 + /3) - 1/3 (2k + 4 + (3) 2 
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